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We consider a diffusive S-N-S junction with electrons in the normal layer driven out of equilibrium 
by external bias. We show that, the non-equilibrium fluctuations of the electron density in the 
normal layer cause the fluctuations of the phase of the order parameter in the S-layers. As a result, 
the magnitude of the Josephson current in the non-equilibrium junction is significantly supressed 
relative to its mean field value. 
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When an electric current flows through a metalhc sam- 
ple, electrons in the metal can no longer be considered 
as a system in equilibrium. In particular the electronic 
distribution function /(e) differs considerably from the 
equilibrium Fermi distribution. The low tenrperature 
two-step distribution function was recently observed in 
tunneling experiments [|lj. 

The ability to change the electronic distribution by 
simply applying voltage to a metallic system opens a 
new possibility to control the supercurrent in a S-N-S 
junction g. The idea is based on the description of the 
supercurrent flow in terms of the electronic states and 
their occupation probabilities /(e). As a function of the 
phase difference of the two superconducting layers in a 
S-N-S junction 9 = 9i — O2 the supercurrent density can 
be written as |3| 



(1) 



J{e)= /de/(e)j(e,( 



where j{e,9) is the contribution of the states with en- 
ergies between e and e -|- de. Changing the electronic 
distribution, one changes the probabilities /(e) and thus 
can change the magnitude g and even the direction Q of 
the supercurrent (in other words, making a 7r-junction). 
Some particular cases of j(e, 6*) were considered theoret- 
ically g. 

It is important to realize that Eq. (p was obtained 
within the mean field approximation where the phase 9 
is fixed. Quantum fluctuations of the phase of the order 
parameter do not induce dramatic change in the ampli- 
tude of the Josephson energy, if the conductance of the 
normal layer is large. The purpose of this Letter is to 
show that it may not be the case in the non-equilibrium 
situation. 

We show that the non-equilibrium fluctuations in the 
normal layer cause the fluctuations of the phase of the or- 
der parameter in both S-layers, thus affecting the Joseph- 
son current in the junction [the term j(e, 9) in Eq. (1^)]. In 
particular, the critical current (which is proportional to 
the Josephson energy Ej) becomes strongly suppressed 
by the non-equilibrium. This effect accompanies the 
change of the magnitude and even of the direction of the 
supercurrent, which is due to the term /(e) in Eq. (Q). 



Thus the critical current in the 7r-junction is smaller than 
the mean field value, as was observed in the experiment 
y|. Here, we will present the phenomenological deriva- 
tion which yields the same results as a calculation based 
on the Keldysh technique (for similar consideration of 
S-N junction see Ref. §]). 

To describe the non-equilibrium fluctuations we first 
identify the collective modes in the junction. For sim- 
plicity we will consider an S-N-S sandwich, where each 
layer can be considered as a 2D film. The effects of the 
finite thickness of the layers will be discussed later in the 
Letter. 

Consider a superconducting filnr at zero temperature. 
All of the excitations with the energy smaller than the 
superconducting gap A are associated with the phase 9 
of the order parameter M. In the isolated S-layer, the 
longitudinal phase fiuctuations correspond to the usual 
2D plasmon with dispersion lo ~ VQ- However, when 
a layer of normal metal is present, the collective mode 
with the linear dispersion relation appears [pj. In the 
S-N-S junction there two such modes corresponding to 
each S-layer. 

The time evolution of phase 9 is governed by hydrody- 
namic equations, which in the absence of external mag- 
netic fields can be written as 
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(2a) 



(2b) 



(2c) 
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where the superscript i = 1, 2 labels the layer, rig is the 
perturbation of the carrier density in S-layers, js is the 
supercurrent, and (p is the electric potential. We wrote 
the London equation (2b) for a dirty superconductor and 
expressed the superfiuid density through the diffusion co- 
efficient Ds in the normal state of the superconductor 
and the thermodynamic density of states per unit area 
in the superconductor Vg . Also in Eq. (0) we neglected 



the terms which describe the Josephson current. We will 
discuss this point later. 

The electron density is normal metal is governed by 
the continuity equation and the Ohm's law 

eiiyn + V • jm == 0; jm = - AnV {e^Vynip + erim) , 

(2d) 

where «„ is the carrier density, jm is the current and Dm 
is the diffusion coefficient in the N-layer. 

Fluctuations of the densities in three layers are coupled 
through the Coulomb potential 

^= fdr'V{r~r')\ni^\r')+n^^\r')+n,nir')\; (2e) 



V = 



(2f) 



So far we have neglected the thickness of the metallic 
layer wherefore only the sum of the electron densities in 
two S-layers is coupled to the density of electrons in the 
metal. Therefore in a strictly two-dimensional model of 
the S-N-S junction with two identical S-layers the Joseph- 
son current (which depends on the phase it difference) is 
not affected by the density fluctuations in the metal. To 
couple the fluctuations in the normal metal to the phase 
difference one needs to introduce some asymmetry into 
the model either by having two different S-layers or by 
taking into account the non-zero thickness of the metallic 
layer. Here we chose the former (see Eq. (pq)). The final 
results do not depend on the asymmetry explicitly, and 
thus are independent of this choice. 

The requirement of the consistency of Eqs. (||) gives 
two acoustic branches of the collective mode correspond- 
ing to the sum and difference of the electron densi- 
ties of the S-layers with dispersion relations wi,2(Q) ~ 
uj'2 — i^i2- The lifetime of both modes is finite. These 
modes are similar to the Schmid-Schon mode in the vicin- 
ity of the critical temperature |g] . The only difference is 
that the normal excitations are not thermally activated 
in the superconductor itself but rather exist in the nor- 
mal metallic layer close to the superconductor, however, 
it does not change the charge dynamics. In the latter 
(odd) mode finite lifetime appears only due to the asym- 



metry of the sandwich 5D = [DT' - DT']/2 



(1) 



1(2)1 



Q 



ttAD, 



'^m{6Df 



A 

¥' 



(3) 



where v^ 



,(1) 



,(2) 



and G, „> denote dimension- 



less conductances of the superconducting (in the nor- 
mal state) and normal layers respectively: Gs,m = 
2TThas,m/ ^"^ = '^T^^i^s.m.Ds^m- The couductanccs are mea- 
sured in units of e^/2nh = 1/(25. 8ifil). The even mode 



is similar to the "phason" mode found in the NS junction 
g . The lifetime of this mode is due to the coupling with 
the relaxation mode in the N-layer 
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Equations (H) and (H) are valid for a;,' > w". For the 
phason mode this condition is satisfied already at small 
frequencies Tllu ~ fiL02 ~ A{Gs/Gm) ^ A. The condition 
for the odd mode is weaker, given the smallness of the 
asymmetry parameter (hvmSD)^ ^ GgGm- 

Now let us consider what happens, when a dc- current 
is driven in the normal layer. The average currents in the 
metal are accompanied by the fluctuations known as the 
shot noise. Since the currents in the metal are coupled to 
those in the S-layers, it is natural to expect that in the 
superconductors the fluctuating currents appear as well, 
and consequently, the phasons are generated. 

To include these fluctuations in our description of the 
S-N-S sandwich, we add Langevin sources 6ji to the cur- 
rent in the normal metal. Equation (Eq) takes the form 
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The Gaussian fluctuations Sji are described by their cor- 
relator. Out of equilibrium, when the energy relaxation 
is negligible t^ —^ 00, the electronic distribution function 
/(e) in the normal metal is the two-step function 
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where ri{x) is the Heaviside function. In that case the 
non-equilibrium part of the correlator of the fluctuations 
can be written as 101 



.•/3 



1 



{Sjt'5j'i)^.Q = -6ai3e D.mVrn{eU - h\uj\)ri{eU - fi\uj\). 

(7) 

The difference of superconducting phases of the S- 
layers 9 = 9i — 62 va. the presence of the current fluctua- 
tions Sji can be determined from the system of Eqs. (||), 
and (0) in the first order in asymmetry 
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Therefore, the correlator of the phase fluctuations has 
two well pronounced poles corresponding to the two col- 
lective modes in the S-N-S sandwich and is proportional 
to the applied voltage U. 

To find the effect of the phase fiuctuations on the 
Josephson current, we need the phase fiuctuations in a 
single point. With the help of Eqs. (M) and (H), we find 
in the lowest order order in asymmetry 
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(9) 

at ?i|a;| < eU and {SO"^),^ = at h\uj\ > ell. Equation (||) 
is valid, provided Ucu » Gg^/Gm- 

Note that in Eq. (0) the asymmetry parameter SD have 
disappeared and therefore only the subleading terms de- 
pend on the asymmetry parameter SD. The leading term 
Eq. (^) is the contribution of the odd mode Eq. (^) to the 
integral in Eq. (Bh. This term dominates because when 
both modes Eq. (^) and Eq. (|j) are well defined, the life- 
time of the odd mode is always longer than the lifetime of 
the even phason mode (provided the asymmetry parame- 
ter 5D <C Dm, Ds). However, the asymmetry parameter 
can not be arbitrary small, as the lifetime of the odd 
mode should be less than the escape time (which is the 
time it takes to remove a phason from the system and 
thus is the characteristic relaxation time). As we con- 
sider the S-N-S sandwich of the infinite size, the escape 
time is practically set to infinity, which is the reason the 
asymmetry disappears from the odd mode contribution 
Eq. (^). In this case, the correlator Eq. (^ is independent 
of the particular choice of asymmetry. The difference be- 
tween various asymmetry realizations becomes important 
if one considers a situation when the odd mode is still 
ballistic, while the even phason mode is already damped. 
We will not discuss this situation in this Letter. 

We have found that in the presence of the normal layer 
the phase fluctuations in the superconductor are large 
due to the large number (~ eU/hio) of the phasons. Now 
we are interested in the effect of these fluctuations on the 
Josephson current. In equilibrium it is determined by the 
difference of time-independent phases of the supercon- 
ducting order parameter in two S-layers J{9) ex sin(6''-*''). 
In the non-equilibrium situation one has to average over 
the phase fluctuations so the Josephson current becomes 
modified by a phase factor. 

The supercurrent flow in Josephson junctions was an- 
alyzed by many authors MM- Equation (nf) expresses 
the Josephson current density as a sum of contributions 
of individual electronic states weighted with their distri- 
bution function. The exact form of these contributions 
j{e,9) depends on the boundary conditions at the S-N 
interface llfl]. In the case of diffusive junction the su- 
percurrent density in the lowest non-vanishing order of 
transparency can be written as lO] 



J(e) = |i?f%sin(0)(e^^«(«)),, 
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(10a) 



(10b) 



where the transverse Thouless energy is Et = fiDn/cf {d 
is the width of the normal layer) and the overall scale is 



given by the bare Josephson energy Ej ' ~ GiG2/vm, 
where Gi and G2 are the tunneling conductances (in 
units of e^ 1%) per unit area of the two N-S interfaces 
in the junction. 

Strictly speaking, Eq. ( |lOa| ) is exact only for homo- 
geneous in space phase fluctuation. The effect of the 
inhomogeneity can be estimated as hDsQ^ / A, while the 
main contribution comes from the odd phason mode with 
DsQ^ ^ ?iw^/A, and so the correction is of the or- 
der h w^/A^ and is small since we consider frequencies 
huj < eU <t: A. 

In equilibrium /(e) is given by the Fermi distribution. 
There are no divergent phase fluctuations, so the factor 
^gi(59(o)^^ in Eq. ( |lOa| ) is just a number which can be in- 
corporated in the definition of Ej . At low temperatures 



the dominant contribution to the integral in Eq. ( |10b[ ) is 
due to small frequencies. At Et ^ T the critical current 
is given by 
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(11) 



In the non-equilibrium situation the critical cur- 
rent ( pT| ) is modified by two effects. First, the distri- 
bution function deviates from the Fermi distribution. At 
low temperatures it is given by the two-step function m). 
The applied voltage eU sets the lower limit for the inte- 
gration in Eq. ( |10b[ ) and for the real part of the integral 
we obtain 
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(12) 



Increasing the applied voltage changes the sign of the 
logarithm in Eq. (p^), see dashed line on Fig. |l|, corre- 
sponding to the TT-junction. 

The second effect is the appearance of the fluctuation 



phase factor in Eq. (10a). To evaluate the average over 9 
we use the phase correlator Eq. iM), Fourier transformed 
to the time domain 
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(13) 



The frequency integration in Eq. ( |13D diverges logarith- 
mically {S9'^) ^ ln(eC//e*). At large frequencies it is cut 
off at TiLu — ell, because at w > ell there are no classical 
fluctuations, see Eq. (O). 

The infrared divergency in Eq. dl3) is due to the fact 
that we have neglected the Josephson current in the con- 
tinuity equation (Ea). When taken into account, it leads 
to the opening of the gap in the spectrum of the collec- 
tive modes, e* = ^/Ej/vm (where Ej is the Josephson 
energy per unit area). This gap, provides the infrared 
cut off in Eq. (O) . It is essential that in the non equilib- 
rium situation, the Josephson energy decreases with the 



increase of the bias voltage U and vanishes at the critical 
point. The actual value of Ej at a given U should be 
determined from the self-consistency equation [obtained 
by substitution of Eq. (O) into Eq. ([TOa|)] 



(14) 



Solving Eq. (|14|) we obtain 




Ej{U) = EfL 
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(15) 



where a = eU jlixt^Gs. The expression Eq. ( [15|) is valid 
when eU > y^Ej/vm- In this case the resulting Ej be- 
comes suppressed by the non-equilibrium (compared to 
its mean field value Ej Lj). 

We now can write down the non equilibrium critical 
current as 

2e 



J. 



-EjiU), 



(16) 



where the renormalized Josephson energy is given by 
Eq. (p^). The dependence of the critical current on the 
bias voltage is illustrated on Fig. ^ 
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FIG. 1. Plot of the Josephson energy Eq. ( |l5| ) as a func- 
tion of the applied voltage U. The dashed line is the volt- 
age-dependent logarithm Eq. (n2|) , which represents the criti- 
cal current without taking into account the phase fluctuations. 

The expression for the Josephson energy Eq. ( [l5|) is the 
main qualitative result of the paper. It shows the effects 
of non-equilibrium are not limited to the change of the 
electronic distribution function. In addition, one has to 
take into account the fluctuations of the superconducting 
phases in both S-layers of the S-N-S junction. The phase 
fluctuations result in two observable effects. First, the 
Josephson energy at the critical point is not described by 
the mean field power law, but exhibits the non-analytic 
behavior Eq. (Eq), illustrated by the inset on Fig. H. Sec- 
ond, when the bias exceeds the critical value (so that the 
critical current becomes negative) the Josephson energy 
is further suppressed relative to the mean field value. 

We should warn the reader that Eqs. (E3), (fq) were 
obtained for the simplest model of the S-N-S junction, 



namely the 2D sandwich. The spectrum of collective 
modes is sensitive to the geometry of the system, there- 
fore, our results are not expected to describe the exper- 
imental data {e.g. of Ref. Q) in detail. However, we 
have presented a strong evidence that the acoustic col- 
lective modes, which are present in the junction in the 
non-equilibrium, can be observed by measuring the sup- 
pression of the Josephson energy. 

In conclusion, we showed that the non-equilibrium fluc- 
tuations of the superconducting phases in the S-N-S junc- 
tion lead to the non-analytic behavior of the Josephson 
energy at the critical point and to its suppression in the 
region of the negative critical current, providing a possi- 
bility to observe the acoustic collective modes (phasons). 

We acknowledge helpful conversations with L.I. Glaz- 
man. LA. is A. P. Sloan and Packard research fellow. 
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